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ABSTRACT 

We investigate the neutron stars spin evolution (breaking, inclination angle evolution 
and radiative precession), taking into account the superfluidity of the neutrons in 
the star core. The neutron star is treated as a two-component system consisting of a 
"charged" component (including the crust and the core protons, electrons and normal 
neutrons) and a core superfluid neutron component. The components arc supposed 
to interact through the mutual friction force. We assume that the "charged" compo- 
nent rotates rigidly. The neutron superfluid velocity field is calculated directly from 
linearized hydrodynamical equations. It is shown that the superfluid core accelerates 
the evolution of inclinaton angle and makes all pulsars evolve to either orthogonal 
or coaxial state. However, rapid evolution seems to contradict the observation data. 
Obtained results together with the observations may allow to examine the superfluid 
models. 
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1 INTRODUCTION 

To date there is no theory explaining the pulsars emission 
mechanism but it is clear that the basic parameters of pul- 
sars signals such as the width of tire pulses profiles depend 
on the spatial orientation of at least two vectors: the angular 
velocity vector and the magnetic moment vector m which 
determines the emission direction. During the neutron star 
life the magnitudes as well as the relative orientation of this 
vectors can evolve. 

Pulsars rota tion dynamics is a c omple x process. As it 
was pointed by iGoldreich fc JulianI (|l969l '). neutron stars 
must possess a very large magnetosphere. Strictly speaking, 
the torque acting on the neutron star can be calculated only 
by using the self-consistent theory o f the magnetosph ere 
which is far from complete at present (|Spitkovskvll2008l ) . It 
makes researchers use the various model expressions for the 
neutron star angular momentum loses torque. The simplest 
one is the radiation-reaction torque, act ing on the ma gne- 
tised sphere which rotates in a vacuum (|Deutschlll955l h In 
that case, the combination P/cosx, where P is the pulsar 
rotation period and x is the inclination angle (the angle 
between the spin axis and the ma gnetic moment), remains 
constant during the pulsa r lifetime (|Davis fc Goldsteinll 19701 : 
iMichel fc Goldwirelll970l ). It leads to a very fast magnetic 
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axis alignment, which appears to contradict the observa- 
tions. This result holds for a fluid star whose rotation defor- 
mation follows the instantaneous rotation axis. 

However, this contradiction seems to be a result of the 
excessive idealization. If the neutron star shape is not per- 
fectly spherical and it is not symmetrical about the magnetic 
axis, the precession c aused by star o blateness will prevent 
the pulsar alignment (|Goldreichlll970l ). 

Except the magnetic-dipolar angular momentum loses 
mechanism pulsars have another one related with the charge 
particles currents in the open field lin es regions ( pulsar 
tubes) of the neutron star magnetosphere jjonesll 19761 ). Cur- 
rents mechanism forces the pulsar to counteralign, and if it 
is comparable with the magnetic-dipolar mechanism, ow- 
ing to their competition the inclination angle evolves much 
more slowly than it fo llows from the vacuum approximation 
iBarsukov et al.ll2013l ). 

In the presence of a magnetosphere the radiation- 
reaction torque itself may, in principle, differ significantly 
from the vacuum case. For example, in t h e fram ework of 
the force-free approximation iBeskin et al.l (| 19831 ) have ar- 
gued that if the tubes currents are absent, a pulsar does not 
lose its angular momentum through the magnetic-dipolar 
mechanism at all. 

Neutron stars are not rigid perfectly, so the different 
dissipation mechanisms may take place in their interiors. 
Because of the finite rigidity the star crust changes its 
shape during the precession period. It leads to the damping 



© 0000 RAS 



2 D. P. Barsukov, 0. A. Goglichidze and A. I. Tsygan 



the precession motion an d decreasing the precession angle 
IChau fc Henriksenlll97ll : lMacv.,1974 ) . 

Neutron stars cores consist mostly of the neutrons with 
small fractions of protons, electrons and probably more ex- 
otic particles like hyperons, kaons, etc. The neutrons in some 
regions of the star should be in superfluid state. Superfluid- 
ity may appear in the rotational dynamics in several ways 
at several time-scales. First, superfluidity is believed to be 
responsible for pulsars glitches, the sudden changes of the 
observational pulsars periods with subsequent smooth recov- 
ery to the initial value. The recovery time-scal es are in the 
range o n days to months. As it was proposed bv lBavm et al.l 
l|l969al ). the central role here plays the pinning and sudden 
unpinning of the vortices governing the superfluid rotation 
to the neutron star crust nuclear lattice. 

Th e vortices pinin g may also cause a neutron star pre- 
cession. IShahamI (|l977l ) has shown that for the neutron stars 
with low oblateness the precession frequency is given by 
ojp = ^'f2, where ^ is the ratio of moment of inertia of the 
pinned superfluid to the rigid component of the rest of the 
star. The such precession must be very fast. The observa- 
tional values of the precession periods obtained for several 
pulsar, contrary to the theory, all are of the order of a year. 
The more detailed theories including the u nperfect pin ning 
do not make the situation more clear (see iLinkI (|2006h for 
the detail). 

Long-term rotation dynamics taking into account the 
internal structure of the neutron stars was studied by re- 
searche r s too. Crust-core coupling has been investigated by 
lEassonl (|l979l ). who has shown that the crust and the core 
plasma are very close to co-rotation which can be supported 
either by the internal magnetic field or by the viscosity. It 
should be noted, however, that he neglected the interaction 
between the core proton-electron plasma and the superfluid 
neutrons. 

The influence of core neutrons superfluidity on th e 
spin evolution was studied bv lCasini fc Montemavorl (|l998l ). 
They treat a neutron star as a two component system the 
components of which rotate rigidly with angular velocities 
Qc and Oq. Crust component interacts with core by the 
friction-like torque iV — — r^^(f2c — fio) and feels the action 
of an external torque. Th e similar problem was investigated 
by iBarsukov et al.l (|2013l ) but they used the more correct 
expression for A'" obta ined directly by the in tegration of the 
mutual friction force (jSedrakian et al.lll999l ). 

The present paper considers the development of the 
semi-hydrodynamical formalism (only superfluid neutrons 
differential rotation effects are taken into account) describ- 
ing the spin evolution of neutron stars which allows to in- 
clude in the equations the realistic mass densities and in- 
teraction coefficients profiles a s well as the different exter- 
nal torques. This work exten ds ISarsukov et al.l (|2013h (and 
ICasini fc Montemavorl (| 19981 )) and it is a step toward the 
full hydrodynamical formalism. 

The paper is organised as follows. In section [2] we re- 
view the basic features of the superfiuid mixture hydrody- 
namics. The proposed semi-hydrodynamical model is de- 
scribed in section [3| Here, we obtain the expression for non- 
axisymmetric core neutrons superfluid flow. The equations 
of motion (described the star breaking, inclination angle evo- 
lution and radiative precession) are formulated in section [l] 
In section[5]we apply some particular models of the neutron 



superfluidity, neutron star core constituents interaction and 
the particular model of the external torque. Obtained results 
is discussed in section [S] 



2 SUPERFLUID HYDRODYNAMICS 

It is well known that for the uncharged superfluids the ve- 
locity held V must satisfy the equation rot v = 0. How- 
ever, the superfluids can rotate by forming an array of vor- 
tices. The velocity he ld near the each vortex has the form 
l|Tillev fc TiUevll 19901 ): 



where mjv is the mass of nucleon forming the superfluid, f„ 
is the distance from the vortex core, ej,^ is the azimuthal 
unit vector. Vorticity is contained only in the vortex cores, 
in which the superfluidity breaks: 

mtv = J2 I ^^kS (fl,, (0 -r)dl (2) 

(|Soninlll987l ). Here, Ck is the unit vector pointing vortices 
orientation, r^^ (l) is the i-th vortex line equation and S{r) 
is the three-dimensional Dirac delta function. Each vortex 
carries the quantum of circulation x = 2irh/2mN ■ 

However, if the characteristic scales much larger than 
the intervortex space Z„ (long wa velength limit), an averag - 
ing procedure can be applied (see lBavm fc Chandler] (|l983l ) 
for the detail). One can introduce a smooth vorticity field 
LOs = n„>fefe, where n„ is the number of vortices per unit 
area. The conservation law for ujs has the form 

dti^s + rotps X Vl] = 0. (3) 

Here, vi is the velocity of the vortices. One can formally 
introduce a vector field Vs such that 

Cjs = rot Vs . (4) 

Substituting this into eq. Q, it is easy to verify that the 
last is equivalent to the Euler-like equation 

dtVs -\- [vs ■ V)vs -\- V/i = fv, (5) 

where 

fv = -UJs X {Vl - Vs) (6) 

and fl is some scalar function. Note that Vs is given up to 
gradient of an arbitrary function. The replacement Vs 
■Us + Vi? leads to replacement /x — >■ /I — ft??— ■ Vi?— (Vi?)^/2. 
Equation ([2} allows us to uniquely determine Vs as < v >, 
where <> denotes the averaging over unit area. In this case, 
function /2 should be a chemical potential per unit mass 
(|Bavm fc Chandleij|l983h . 

In order to exclude vi from the equation let us consider 
the force balance for a vortex. On the one hand, if the super- 
fiuid coexists with some liquid which moves with velocity Vc 
(index "c" is used in order to unify the notations with the 
following sections of the paper), the vortex feels the action 
of the friction force. A lot of friction mechanisms can take 
place in the neutron star interior but, in any case, if the 
friction is not very strong, the force acting per unit length 
can be represented as 

F,^~^{vi-v,), (7) 
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where pc is density of "c"-liquid, Tc is the characteristic time- 
scale depending on particular friction mechanism. If there 
are several comparable mechanisms in some star region, Tc is 
calculated by the formula = '''ci +'''c2 ~^ ■■■'^'''ck ^ where 
Tci are the time-scales specifying each mechanism separately. 

On the other hand, if the vortex moves through the 
superfluid, it feels the action of the Magnus force which has 
the form 



Fm, = pa>i:ek X {vi - V - rot A). 



(8) 



Here, rot A is the term arising from th e vortex local self- 
acting (see, for example. ISchwar"3 (|l985l ')'l. A is a vector par- 
allel to el. 

The thickness of the vortex core equals to the coherent 
length 5„ which for the ne utron star core is o f the order 
of several tens nucleon sizes (|Bavm et al.lll969bl '). The mass 
contained in the vortex core are negligibly small. Thus, New- 
ton second law for the vortex reduces just to 



Fm Fc 



(9) 



Solving equation Q for vi, averaging the result over 
unit area and substi t uting it into Q, o ne can obtain 
l|Sedrakian et al.lll999l : [Siderv fc AlpailEoogl ') 



/„ = —uja X rot A„ + P'uis X {Vs - 
+I3ek X [ids X {vs - 

where we introduce two coefficients 



/?' = 



Vc + rot A„)-|- 
- + TOt X„)] , (10) 



(11) 



where 



Pc 



pa TcUJa 



(12) 



is the coupling parameter. This kind of interaction taking 
place in superfluids is usually called the mutual friction 
iHall fc Vinenlll956l ). Depending on a there are two regimes. 
Weak-coupling regime corresponds to a ^ 1. In this case, 
~ P' 1 and vi ~ Va- If the opposite inequality a ^ 1 
satisfies, the strong-coupling regime takes place. It means 
that /3 ^ 1, /3' « 1 and vi « Vc- Note that this analysis 
is quite general. The only thing which determines by the 
physical nature of "c" -liquid is Tc- 

So, in the long wavelength limit there is no need to con- 
sider the vortices dynamics. Instead that one can solve the 
Euler-like equation to which the mass, momentum and en- 
ergy conservation laws for the whole liquid (including ther- 
mal excitations) should be added of course. 

The hydrodynamical equations for the rotational su- 
perfiuid liquid s can be obtained f rom the conservation laws 
consideration (|Khalatnikovll2000l '). This method don't deal 
with the vortex lattice at all, but it is quite rigorous and 
it naturally can be extended to the more general superfluid 
systems. One just need to assume that the internal energy 
depends on the Cj. 

An example of such systems i s a mixture of superfluid 
liquids. lAndreev fc BashkinI (|l975l ) have shown that the su- 
perfluid hydrodynamics must be modified in that case. The 
strong interaction between the particles forming the two su- 
perfluid constituents gives rise to the so-called entrainment 



effect which leads to the following relations: 

Jn — (Pn pnn pnp )Vex + PnnVn + pnpVp, (13) 

Jp — (Pp Ppp Ppn)'^^ex PppVp PpnVnt (1^) 
Pn — (Pn — Pnn — pnp)Vex + {pnn + pnp)Vn, (15) 
Pp — {pp — Ppp — Ppn)Vex + {ppp + Ppn)Vp, (16) 

where Ja is the mass current of the a constituent (a — p,n), 
Pa are the momentum densities, Vex is the thermal excita- 
tion velocity (it is assumed that the excitation velocity is 
the same for both fluids), p^a is the mass density matrix 
(also called the entrainment matrix or the Andreev-Bashkin 
matrix). It can be shown that pnp ~ Ppn- So, in the super- 
fluid mixtures the mass currents corresponding to each fluid 
are no longer parallel to they velocities and don't equal to 
corresponding momentum density. Note, however, that 



Jn ~\~ Jp — Pn H~ J 



(17) 



The superfluid velocities here are determined just like for 
the ordinary superfluids: 

Vn = -^VSn, Vp = TT^VS'p —A, (18) 

2m„ Inip nipC 

where Sa are the phases of the complex order parameters, 
A is the electromagnetic vector potential. 

One can also introduce the mass currents and momen- 
tum densities for the superfuids and thermal excitations sep- 
arately: 



Ja = 

Ps = psVn, 
Ppis) = Pp{s)^p, 
Pex — pex'^e.x-) 



(19) 
(20) 
(21) 
(22) 
(23) 



where we denote pa = pnn + pnp, Pp{a) ~ Ppp + Ppn and 
pex = pn + pp ^ Pa — Pp{a) ■ Here aud after we will use the 
index "s" instead "7i(s)" for all quantities relating to the 
superfluid part of the neutron liquid. 

The system of hydrodynamical equations describ- 
ing the mixture of superfluid neutrons, superconduct- 
ing protons and degenerate electrons and muons taking 
into account the neutron-proton entrainment and gravi- 
tational force was develop e d by Mendell and Lindblom 
(|Mendell fc LindblomI [l99ll : iMendelll [l991al lbh. They used 
the Khalatnikov's method. The "superfluid neutron" part 
of Mendell-Lindblom system of equations has the form 



dtVa + {Va ■V)Va+V{p + $G) 



pnp 



{Vp~Va)Xi^a+fv, (24) 



dtpa + div Ja = Fa 



fv = iOa X rot A„-|- 

Pa 



(25) 

(26) 



+ 13' 
+ pSk X 



I Pnn , pnp -, ^ , 1 ^ r* 

LJa X Va H Vp — «e H rOt A„ 

Pa pa Pa 



, Pnn -* , pnp -* , 1 , T* 

iJa X I Va H Vp — Ve -\ rOt A„ 

Pa Pa Pa 



Here, r„ is the superfluid neutrons creation rate. An ~ 
dUo/duJa, Uo is the neutron liquid internal energy density. 
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Ve is the electrons velocity. It should be noted that the phe- 
nomenological approach by itself doesn't allow to estimate 
the mutual friction coefficients. Here, j3 and /3' are just the 
arbitrary parameters. 

Mendell and Lindblom used the Khalatnikov's phe- 
nomenological approach but the same equations can be ob- 
tained from the vortices dynamics consideration if we sup- 
pose that the Magnus force takes the form 

Fra = >i:ek X (^paVi - Js - rot , (27) 

and that the friction force equals to 

Fc = ^{vi-v,). (28) 

The last assumption means that only the electrons scat- 
tering is taken into account. In present paper we restrict 
ourselves to considering only electron-scattering processes 
but in final equations the particular scattering mechanism 

specified only by /3 and coefficients. 

As it follows from the calculations of iMendelll (|l991ah . 

X„^^^\n(^]e,, (29) 

4:71 Ppp VCn/ 

where — p„„ppp — p^p is the determinant of the mass 
density matrix. 




Figure 1. The cylindrical coordinate system and the vectors used 
in section [3] 

It is important to note that by the Uc we will denote the 
quantity .1^1 Pc (as it follows from pop ). so Ua is not equal 
to Up. 

Equations (|24|) - l|26|l in co- rotating frame take the form 

dlus + 2[n X Us] + {us ■ V)us + Vpi = X r\+ (36) 
+ — {up - Us) X (2Q + rot Us) + fv, 

Ps 



3 THE MODEL 

We treat a neutron star core as a two-component system. 
The first component consists of the superfluid neutrons and 
its dynamics is described by equations (|24p -(|26 p . Protons, 
electrons, and normal neutrons are coupled to each other on 
the small time-scales and form the second component which 
we denote as a "charged" component. All particles forming 
the charged component are supposed to move with the same 
velocity Vc- 



Vex — Jp(s)/Pp(a) 



(30) 



The star crust rotates with the angular velocity Q which 
is identified with the observed angular velocity of the pulsar. 
The external torque K acts on the crust. We will consider 
the torques depending only on Q and the configuration of 
magnetic field which supposed to be frozen into the crust. 
In this case, the torque is a very slowly evolving vector 
in the frame co-rotating with the crust, so this frame is the 
most suitable for the considering problem. 

First of all let us introduce the velocity fields and the 
mass currents measured in co-rotating frame: 



ita = t'a — X r], 



(31) 
(32) 

Here, a is the constituent index (a = e,s,p{s),ex). More- 
over, assumption (|30|l allows us to introduce this quantities 
for the whole charged component in a similar form: 



where 



Uc ^ Vc — [il x r\, 

J* = Jc - Pc[^ X f]. 



pc = ^Pa. 



(33) 
(34) 

(35) 



dtPs -l-div J* = Ts, 



fv = (2f2 + rot Us) X < rot A„ + 

nf I Pnn Pnp ^ ^ , 1 j_ \ \ , 

+ p Us H Up - Ue H rot A„ + 

Ps ps ' 



(37) 



(38) 



+/3 



ps 

20. + rot Us 



Pnn 
ps 



Us + 



Pnp ^ 
ps 



Ue H rot A„ 

Ps 



\2Vt + rot Us I 

where dl denotes the time derivative taken in the co-rotating 

frame, jli = -(1/2) [fj x f]^ 4- /i + $g, ^ = dtQ. 

It will be useful further to expand Sl = + Six, where 

^11 = (d-Q.)n/Q.'^ and fix = 0-^11. Based on this expansion 
let us introduce three orthogonal unit vectors Cx, Cy, as 
follows: 



(39) 



Let us also introduce associated cylindrical coordinates f , </>, 
and z {f = Q, z = Q point the centre of the star) with 
corresponding orthogonal unit vectors , and e, . We will 
use the notations: f for the cylindrical radius vector and r 
for the spherical one (see fig.[T]). All densities and coefficients 
are supposed to be functions only of r. 

We assume that the following inequalities are satisfied: 

Slf > Us > Uc (40) 

so we keep in equations (|36p - p8p only the linear in Us terms 
and neglect all terms containing Uc and Ucx ■ For the protons 
the particles (and charge) current is proportional to Jp, so 
we suppose that 



Jp(s)/Pp{s) — PeUc — p. 



ex ^ex 



0. 



(41) 
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Together with ([T9| and ((20l) it leads to 



and 



j: 



Ppp 



(42) 



(43) 



The second part of inequality (|40p requires a sufficiently 
effective physical mechanism damping the differential mo- 
tions of the charged component and its strong coupling with 
the crust, ft can be ensure by viscosity or by the magnetic 
field (the last requires the second type superconductivity 
for the protons). In any case, the more self-consistent model 
should involve hydrodynamical consideration for the both 
components of the star core but it significantly complicates 
the problem. We will include the Vc perturbation in our fu- 
ture developments. 

Let us discuss the rot A„-containing terms of equation 
(|26p more detail. One can write 

dX„ . 



rot A„ 



df 



+ rot[A„(efe — e,)], 



where vector el- — can be expanded in terms of Us 
1 



20. 



X 



[ez X rot Us] + O (ug 



Using (|29|l . one can obtain 

LOs X rot An ~ T^Er 

ps ps or 

where we introduce a parameter 



In 



0{Us€) 



(44) 



(45) 



(46) 



(47) 



'i-ITpsPpp \ Q, J \L J \^n, 

which is much less than unity. Here, L is the characteristic 
scale. In fact, it means that the linear in Us terms of the 
expansion of uja x rot An actually in addition to its contains 
another small parameter e. Hence, for the rot An-containing 
terms we suppose that 



Lds X rot A„ 



or 



(48) 



is a good approximation. From a physical point of view it 
means that we neglect a part of vortices self-acting arising 
from their small curvature. 

Assumption (|48l) is a priori incorrect close to the super- 
fiuidity breaking surface. The vor tices lines should satisfy 
the following boundary condition (|Khalatnikovll2000l ) : 



0, 



(49) 



where Ts is the radius of sphere in which the superfiuidity 
breaks (the radius of the superfiuidity core) and vector n is 
normal vector to this surface. It means that = n on the 
surface and it significantly differs from there. However, 
the departure of ek from z direction exponentially decreases 
depthward the superfiuid core with a characteristic scale 



Ps 



10^ 



10" cm 



(50) 



It easy can be obta ined using eq. (8) by 
ISedrakvan fc Savvi HI 1 19791 ) So, this effect seems to 
be insignificant for our problem. 

Vector Q as well as the external torque K evolves very 



slowly in co-rotating frame in compare with the mutual fric- 
tion time-scales. Thus, the long-term evolution corresponds 
to the quasistationary solution of equations and (f57)) . It 
means that the time-derivative terms d* become negligibly 
small on the large time-scales. 

Continuity equation p7p can be rewritten as 



d'pa + div J* = dtPc + div J*. 



(51) 



Here, 9*Ps and d* Pc equal to zero as we just argued. Charged 
component mass current is negligibly small because of (|4ip . 
Thus, Ts equals to zero up to corrections ~ Uc- 

In present paper we also neglect the gravitational po- 
tential perturbation. 

After all simplifications we obtain the following equa- 
tions for It.: 



20/36*2 X 
-2^7 



^ 1 IdXn^ 
lis --^e^ 

PppPs Ps r or 

^ 1 19A„^ 
Us -^T^e^ 

Pppps Ps r Or 



= [fi X r] 



(52) 



div 



Its 

Ppp 



= 0, 



(53) 



where 7 = (1 — /3'). In addition, the velocity field must satisfy 
the boundary condition 



[Jn ■ r|r 



Ppp 



-(us -r) 



0, 



(54) 



Note that all mass densities in (|52p - (|54|l are not perturbed. 
The solution satisfying (|52p , (|53p , and (|54p has the form 



.BizP±[axr\ 



pppps 7 



+ e. 



ppp ^11 
f|2 20 



1 d 



20. 72 + 



(55) 



+ 



»|| ;fl + 7V' 
20 72 + /32 

7 

72 4-/32 



fdf V ^'72-^/32 



/3 fix 



dz' , 



72+^2 Q 



(56) 



O I 



, 1 19A„^ 

H --^e^ 

Ps r or 



pi 



yzO± — f^ii / ip{r )f dr , 



:dz' 



:dz' 



^2+^2 



= \/r 



(57) 

(58) 
(59) 



First, let us consider an uniform rotating star {O — 0) 
In this case, the velocity field reduces just to 

ppp dXji 



Us = ^-^e^ 



(60) 



but it is not equal to zero. It means that even in station- 
ary situation the superfluid motion do not relaxes to rigid- 
body rotation. However, there is no paradox here. One just 
needs to recall that the vortices are the things which interact 
with the charged component but not the superfluid itself. It 
means that the friction tends to make vi equal to Ve- If we 
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return to vortex force balance and suppose that vi — Ve, we 
obtain 



Pnn -t Pnp -* 

ps>i:ek X \vi Vs 

Ps 



— rot A„ = 

ps ' Ps J 

---^{vi-v.) = 0, 



Pnn -* Pnp 1 j_\*iT^^ 

Vs H Vp — Ve = rot \n + Mlfifc, 

Ps ps Ps 



(61) 



(62) 



where V\\ is an arbitrary function. Taking into account 
equations (|42p and (|48p . we obtain exactly equation (|60p 
(V^l equals to zero because of continuity equation (|53|l and 
boundary condition ((54}). So, the mass density gradient in 
star core gives rise to the force pushing vortices from the 
rotational axis. 

In general case, the flow has a complex form. It con- 
sists of a differential rotation about the local axis -w and a 
poloidal flow. The last, however, is the consequence of the 
variation of the mutual friction coefficients with r and it 
vanishes if we make /3 and /?' to be constant. If we, in ad- 
dition, neglect the entrainment effect [pppps/o^ ~ 1) and 
the vortices energy (An ~ 0), we obtain a rigidly rotating 
solution. This simple ca se has been studied b y the authors 
in their previous paper (jBarsukov et aLlbOlSl ). 

The inequality Us <g Qr must be satisfied for the solu- 
tion to be valid. If /3 and /3' can be expressed in terms of 
coupling parameter a according to pip , taking into account 
that in this case 



7 



/3 



7^ +/32 



we obtain the following conditions for a: 



(63) 



(64) 



On the one hand, the vortices lattice should interact with 
charged component efhciently enough to rapidly reform with 
respect to the changing of Q. One the other hand, if the 
interaction will be very strong, the vortices can not move 
through the charged component almost at all. The absence 
of the interaction as well as the very strong coupling leads 
to that the difference between Vs and Vc during the rotation 
evolution can become not small. 

Protons of course are not required to be supercon- 
ducting everywhere inside the superfluid core. Entrain- 
ment effect vanishes whe re the superconductivity breaks 
(IGusakov fc Haenselllioosi ). Obtained solution, however, re- 
mains valid there. One just needs to replace 



PppPl^l and 

Ps 



(65) 



Except this formally replacements the superconductivity 
breaking apparently leads to the mutual friction coefficients 
decrease by several orders of magnitude (see the detail be- 
low). 

The solution formally tends to infinity on the superfluid 
breaking surface. It means that Ma can not be treated as a 
small perturbation in the boundary region and one needs to 
solve the non-linear hydrodynamical equations there. How- 
ever, the thickness of the layer in which the non-linear terms 
become important seems to be very small. Moreover, the su- 
perfluid mass current which we will use in next section to 



formulate the equations of pulsar motion obtained with (|55p - 
H59p tends to zero there. So, this layer seemingly does not 
play a signiflcant role in the angular momentum transfer. 



4 THE EQUATIONS OF MOTION 

Possessing the expression for the superfluid velocity field, 
one can calculate the angular momentum transfer rate from 
superfiuid core to the charged component. From the angular 
momentum conservation law we have 



M = 



(66) 



where M is the full angular momentum of the star, which 
obviously is just a sum of angular momenta of the compo- 
nents. So one can write 



fx (pc +Ps)dV, 



(67) 



(68) 



The taking into account H4ip and (|43p after some integrating 
leads to 



M = Af, + Ms ■- 
and, using (fT7)). 

M ^ f fx {Jc + Js)dV. 



-Is^ 



f Ppp 



dus 
dt 



dV + K, 



(69) 



where we introduce the moments of inertia for the compo- 
nents 



stt r' 

3 J 





where ps and pc are the mass densities of the superfluid and 
charged components, Icr is the crust moment of inertia. 

The next step is to differentiate the velocity field Us 
with respect to time. Doing this, we as before neglect dtUs 
term, so 



dtUs ^ [Q X Us] ~ {[O. X r\ ■ V 



(71) 



Substituting ([55}, ([56} into ([TTJ, and ([TTJ into ([691), we ob- 
tain 



Jef^ = Sllsfl - Sils^W + S3ls[e, X + K, 

where we introduce the following coefficients: 

Stt ? 7/?' - fi^ 



Si 



3Is J 7^ + 





■psT^dr, 



7 



psT^dr, 



8tv_ 

Sis J 7' + /?' 



^ psr'^dr, 



(72) 

(73) 
(74) 
(75) 



Note that if expressions pip are correct. Si equals to zero 
and S2 equals to unity. Equation S2 — Si = 1, however, is 
satisfied without any assumptions about the relations be- 
tween /3 and /?'. 

Next, we need to introduce a new basis e^, Sy, and e^, 
where Ez — rh/m and vectors Ex and Ey are anchored in the 
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star crust and perpendicular to and to each other. In this 
basis the orientation of Sl determines by two angles x s-nd 
ipn (see fig. [2]). Angle x is the pulsar inclination angle, and 
the variation of (fin relates with star precession. Without 
making any additional restriction external torque K can be 
represented as 



K = Kne^ + K^t, + Kj_[e, x 



(76) 



Equation (|72l) can be solved for Ci and rewritten as three 
scalar equations which describe pulsar braking, inclination 
angle evolution and torque-driven precession respectively: 

Kn + Km cos X 



n = 



Is+h 



X o 



n (Je-Sl/.)2+. 



51 Jf 



smx, 



1 (Jc ~ SiL)Ka_ + S3lsK„ 

n (/, - + 5|72 



(77) 



(78) 



(79) 



Equation (|77[l formally has the sam e form as the one for 
the rigid star (jBarsukov et al.ll2009l ). Neglecting d*Us, we 
assume, in fact, that the energy loss of the superfiuid com- 
ponent equals just to {Is/2){Q-Cl). However, it doesn't mean 
that the breaking will occur with the same rate. Equations 
(|77p - (|79p should be solved simultaneously. The presence of 
the superfluid core leads to that the torque perpendicular 
term K± starts to affect the inclination angle evolution (see 
ICasini fc MontemavoJ (| 19981 )) as well as term Km starts to 
affect the precession. 



5 APPLICATION 

Up to now the developed formalism does not require any 
specification of the profiles of the mutual friction coefficients 
and the mass densities as well as Kn, Km and K± can be 
arbitrary small in magnitude slow varying functions. Next 
we apply the particular models (simplistic at some points) 
in order to demonstrate that the choice of the model may 
significantly affect the rate of the neutron star rotation evo- 
lution. 

Mass density matrix pa/} with the taking into ac- 
count the temperature effects has been calculated by 



iGusakov fc ^Haenselliooi). We do not give here the lengthy 
expressions obtained by Gusakov and Haensel. Note, how- 
ever, that the densities pa/s depend on the parameters = 
T/Tca, where T is the mixture temperature and Tea are the 
critical temperatures of proton and neutron fluids. Gusakov 
and Haensel have shown that the magnitudes of paa and 
Pap tend to zero with T approaching Tea- It means that the 
entrainment effect vanishes when the superfluidity breaks 
for one of the fluids. We have used this in (|65p . 

The critical temperatures themselves depend on the 
mass density. Unfortunately, these dependences calculated 
by different authors m ay significantly differ from each other 
lYakovIev et al.ll 19991 ). It may leads to the large variation of 
the amount of the neutron superfiuid in star cores and the 
mutual friction coefficients. The aim of this section is to in- 
vestigate the sensitivity of the long-term rotation dynamics 
to these variations. 

Despite the differences between the nucleon superflu- 
idity /su perconductivity theories there are some common 
results (jYakovlev et al.l 1 19991 ). Each critical temperature 
proflle should possess a maximum T^"^ at some density 
p^a^ ■ For the protons exceeds the nuclear density 

po ~ 2.8 X 10^* g/cm^. For the neutrons, actually, there are 
two types of the neutron superfluidity. Singlet-pairing (^^o) 
takes place in the neutron stars crusts. We neglect this phe- 
nomenon in present paper. In the star core the microscopic 
theories predict a triplet-pairing type (^^P^) superfluidity for 
which > and TZ^^ < r™»^ 

We model the critical temperature profiles by the 
parabolas (see fig[H|): 



log 



Tea 
'T'maa 



-41og^o 



(80) 



We flx T^"''' to be equal to 5 x 10^ K and variate T™"^. The 
densities p^"'^ and p™"^ corresponding to the maximum of 
the proflles both are put to be equal to 8 x lO'^^g/cm^. It is 
not quite realistic of course but the variation of T^'^^ allows 
us to investigate the sensitivity of the results to the choice 
of the critical temperature proflle. 

We will consider a neutron star with radius r„s ~ 12 
km. Parameter Vs is determined as a radius of the sphere 



on which T = T„ 



{")■ 



The star core is assumed to be 



isothermal. Its temperature does not change with time and 
equals to 5 x 10^ K. In order to calculate pp and p„ proflles 
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we use iHeiselberg fc Hiorth-JensenI (|l999t ) approximation of 
ARP98 equation of state. 

According to (|14|) the neutron vortices must generate 
the magnetic field with typical fluxes which equal to 



PPP / 



$0, 



(81) 



where c&o = 2nhc/2e is the flux quantum (jAlpar et al.ll 1983 ). 
We suppose that the most effective mutual friction mecha- 
nism is based on the electrons scattering on the vortices 
magnetic field. For thi s mechanism the co upling parameter 
has been calculated bv lAlpar et al.l (|f 9841 ') and it equals to 



UALS 



1.3 X 10" 



1 



{xpi4) 



1/6 



\pppj 



(82) 

where x = Pc/p, pu = p/lO^^g/cm^, p ^ pc + ps, m* is the 
effective proton mass. Thus, it is the weak-coupling regime, 
so 



Si =0, & = 1, 53 < 1. 



(83) 



Not that this mechanism requires the protons to be in su- 
perconducting state. In our simplistic model it satisfies ev- 
erywhere inside the superfiuid neutron core. The neutron 
vortices have also the own magnetic field caused by the 
spontaneous magnetizatio n which should ta ke place in 
neutron superfiuid liquids (|Sauls et al.lll983 ). However, this 
field is several orders of magnitude smaller than the field 
generated by proton currents. So, if the superconductivity 
breaks in some core region, the mutual friction become much 
less effectively. The electrons also can scatter o n the ther- 
mal excitation in the vortices cores (fFcibclman| [l97ll ) This 
mechanism is very sensitive to the temperature, however, 
it apparently is much less effective then the Alpar-Lauger- 
Sauls scattering even for the young neutron stars. Except 
the electrons the neutron vortices can interact with the pro- 
ton fiux tubes. However, we will suppose either that the 
superconducting protons are in the first type state or that 
this kind of interaction is much less effective then the Alpar- 
Lauger-Sauls scattering too. 

As for the torqu e K, we use the one proposed by 
iBarsukov etai] (|2009l '): 



- q(x. Vn)) cosxem — + 
+Rcosx[ez X Cm]) 



where 



Ko 



3c3 



R 



9 



10 Vtrn 



(84) 



(85) 



This torque is a sum of magneto-dipolar angular momen- 
tum losses torq ue calculated in the vacuum approximation 
Deutschlll955l ) and the one related with the currents losses 
JoneElll97d r iere, a(x,^n) is the function ~ 1 which de- 



pends on the structure of the small-scale magnetic fields in 
vicinities of neutron star magnetic poles. Small-scale fields 
bend the pulsar tubes and affect the tubes currents. In the 
case of pure dipolar magnetic field, function a(x, V^n) is a 
constant. 

The last te rm in (1841) lea ds to the torque-driven radia- 
tive precession (|Melatoj|2000l ). Note that this term contains 




Figure 4. The behaviour of function ai/(x) for different values 
of non-dipolarity parameter u. 



the big parameter 



i? 3.6 X 10'' 



1 sec 



(86) 



so K_L 'S> Km and Kn for this torque. Together with 
it means that the first term in (|79|l much greater than the 
second and, consequently, 

IcKoR 



■ cos X- 



(87) 



If X is not too close to n/2, the radiative precession 
period 



27r 



Tp = 

KqR Jc cos X 



1.9 X lO'' 



r2 I 02^2 

-"0(45) + '-'3-', 



(45) 



1 sec J Bi2/c(45) cosx 



years 



is three orders of magnitude smaller that the characteristic 
time-scale of the pulsar braking and inclination angle evo- 



lution. Here, L 



(45) 



are the moments of inertia measured in 



units of 10^^ g/cm'^, _Bi2 is the dipolar magnetic field mea- 
sured on the magnetic poles in units of 10^^ G. This fact 
allows us to average the equations over Tp. Dividing aver- 
aged equation (|78p by averaged (I77 |, taking into account 
that P = 2tv/Q,, one can obtain (see IBarsukov et al.l (|2013l ) 
for the detail) 



dx _ 1 Icjlc + Is) sinxcosx 
dP^ P /2 + Sill sin^ X + (x) cos2 X ■ 



(89) 



where a„{x) = {2-k)^^ j'^^^ OL{x,'Pn)d(pa.- The behaviour of 
a( Y) calculated in the framework of the model proposed 
bv IBarsukov et al.1 \2Qli ) is shown in figure (4] In proposed 
model the behaviour of oli,(x) is determined only by non- 
dipolarity parameter u — Bi/Bo, where Bo and Bi are the 
magnitudes of dipolar and small-scale magnetic field mea- 
sured on the neutron star polar caps. 

Equation a(x) = 1 has a root Xeg for a range of non- 
dipolarity parameter v values. For a perfectly rigid star Xeq 
is the equilibrium inclination angle. 

Inclination angle evolution as a function of the rotation 
period does not depend on the magnitude of the magnetic 
field but only on the relation 1/ — Bi/Bq. We suppose that 
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Table 1. The values of Is/{Is + /c) and ^3 for different T™<»^ 



T™"^ (xlO^K) 


Is/(Is+Ic) 


53 (xlO-^) 


1 


0.07 


2.38 


5 


0.33 


5.08 


20 


0.5 


5.05 



non-dipolarity parameter v remains constant during the evo- 
lution. This assumption supported by some magnetic decay 
models according to which the par ameter v evolves much 
slower than the magnetic field itself (|Mitra et al.lll999l l (We 
are interested in the small-scale fields with characteristic 
scale ~ 1 km). However, it should be noted t hat it is not the 
generally accepted result (jUrpin fc Gilll2004h 

The last term of equation (|89p (see also eq. (ITSp ) con- 
tains K± which is not affect the inclination angle evolution 
in the case of rigidly rotating star. The influence of this term 
becomes important when the rotation period approaches to 
the value 



Ic 207rr„ 
'Sjl 9c 



(90) 



So, the greater the amount of the neutron superfluid and 
the stronger the interaction, the sooner it happens. When 
this term becomes dominant, it makes all pulsars evolve to 
the orthogonal state. 

The calculated values of Js/(/s+Jc) and S3 for different 
T^°'^ are given in table [T] The evolution trajectories for 
different initial % ^-nd P and different y™"^ obtained by the 
integration of equation (|89[) are shown in figure [5] 

The observational data on 62 pulsars also are plotted in 
the figure. We use the data from (Malov 20041 ) on pulsars in- 
clination angles measured by using the maximum derivative 
of the position angle (denoted as P2 is the book). 



6 DISCUSSION 

The semi-hydrodynamical consideration results in equation 
(|89p which is equivalent to the analogous equation ob- 
tained under the assumption that the mutual friction co- 
efficients are constant and both components rotate rigidly 
iBarsukov et al.ll2013l ). One just need to replace cr by its 
value averaged over the core volume with the weight function 



^ As one can see from table[T]the value of Si does not 



change significantly with T™"^. However, the higher values 
of y™"^ correspond to higher values of relation Is /{Is + Ic) 
and to the more rapid inclination angle evolution. 

In more general root of equation Km{xeq) ~ 

is the stable inclination angle when ^^^(xeg) > for the 
rigidly rotating star (the bar denotes the averaging over 
ipn)- Even if the rate of the angle evolution is high, the 
observed wide distribution of x can be explained by the ex- 
istence of the equilibrium inclination angles which may be 
different for each pulsar. For example, u should vary from 
one pulsar to another. But now, as it follows from equa- 
tion (|89p . the superfluid core makes the pulsars evolve to 
the orthogonal state and its influence increases during the 
star breaking. So, fast evolution, corresponding to the large 
amount of the superfluid neutrons, seems to be in contra- 
diction with the observational data. This facts may allow to 



0.0 



0.5 



-0.5 

lg(P/l s) 

a 



lg(P/l s) 

b 




★** 



-n.5 



0.0 



lg(P/l s) 
d 

Figure 5. The evolution trajectories obtained (a) for the rigidly 
rotating star {Is = 0) and for the "two components" star with 
T™"^ equal to (b) 1 X 10* K, (c) 5 X 10* K (d) 2 X 10^ K. Dotted 
lines correspond to the initial period Pq equal to 10 msec, and 
dashed lines correspond to Pq = 100 msec. The non-dipolarity 
parameter u is taken to be equal to 0.5. 
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examine the superfluid models. Note, however, that we've 
considered a perfect spherical star. The oblateness along the 
magn etic axis which should be caused by strong magnetic 
field (iGoldreich 1970l) just leads to the redefinition of the R 
parameter ( Barsukov fc Tsyganl I2OI0I '). If the oblateness is 
strong enough, parameter R may become negative, and the 
pulsars will tend to coaxial state instead the orthogonal. But 
it doesn't change the situation qualitatively. Triaxial stars, 
in principle, may possess the equilibrium inclination angles 
different from and 7r/2. This question requires a separate 
studying. 

In section [S] we've supposed that the protons are in 
the superconducting state. If the state is of the first type 
iBucklev et al.ll2004l ). there is no magnetic field inside the 
superfl uid core. In the case of the second type superconduc- 
tivity m et al.l If 969bt ) , magnetic field contains in the 
flux tubes each of which carries the quantum of the magnetic 
flux $0 = hc/2e. The flux tubes density can be estimated 
n/ ^ B/^o = 5 X fO^^-Bi3 cm~^, where B13 is the aver- 
aged over the area ^ ^''"^ magnetic field in units of f O^^G. 
Meanwhile, averaging ((2| over the area 3> l^, one can obtain 
riv ~ 6.4x fO^P~^cm~'^. So, flux tubes much more numerous 
than the neutron vortices. It allows to replace the flux tubes 
lattice by the continuum medium when a vortex dynamics is 
consi dered. If the core magn etic field has the chaotic struc- 
ture (|Ruderman et all 1 19981 ). the interaction with the flux 
tubes can be included in the model just by replacing aALS 
by (JALS + 0' ft, where aft describes the vorte x friction on the 
flux t ubes medium. The calculations show (|Siderv fc AlpaJ 
I2OO9I ) that the interaction with flux tubes much more effec- 
tive {aft S> aALs), so the inclination angle evolution should 
be more sensitive to the relation h/Ic- If the magnetic fleld 
structure has a preferred direction (for exam ple, if it has 
large toroidal component (|Braithwaitdl2009h '). the friction 
becomes anisotropic. That case goes beyond the scope of 
the model proposed in present paper. 

We've assumed that Ms ^ Uc. In fact, it is a strong as- 
sumption, which is not quite established in present paper. 
The more self-consistent consideration requires the inves- 
tigation of the influence of the charged component veloc- 
ity perturbation on the rotational dynamics along with its. 
We do not take into account any pinning phenomena which 
may play significant role in neutron star rotational dynam- 
ics. Also we do not take into account the thermal evolution 
of the neutron stars. We are planing to include this factors 
in our future developments. 
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